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Abstract 

We investigate the production of positive parity baryon resonances in proton electromag- 
netic scattering within the Sakai-Sugimoto modeL The latter is a string model for the 
non-perturbative regime of large QCD. Using holographic techniques we calculate the 
generalized Dirac and Pauli form factors that describe resonance production. We use these 
results to estimate the contribution of resonance production to the proton structure func- 
tions. Interestingly, we find an approximate Callan-Gross relation for the structure functions 
in a regime of intermediate values of the Bjorken variable. 
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A Breit frame for inelastic scattering 



1. Introduction 



In the regime of low momentum transfer (yg^ lower than a few GeVs), non-perturbative 
effects become relevant in hadronic scattering. Lattice QCD (the numerical approach of 
substituting the continuous spacetime by a lattice of points) is quite difficult to use here 
because of the real time dependence in scattering processes. Effective models can be very 
useful but they have a low predictability because they may depend on many parameters. 

In recent years, gauge/string duality has evolved into an invaluable tool to study many 
strongly coupled phenomena in particle and condensed matter physics. In particular, it has 
brought a new insight into the strongly coupled regime of Yang-Mills theories. A good ex- 
ample is the successful prediction of the shear viscosity for strongly coupled theories whose 
gravity dual involves a black hole in Anti-de-Sitter space [1]. This prediction is in accor- 
dance with the small value observed in the quark-gluon plasma phase observed in heavy ion 
collisions at RHIC and LHC. The study of 10-d string models (top-down approach) and 5-d 
effective models (bottom-up approach) dual to 4-d QCD-like theories is known as AdS/QCD 
(see [2l[3lB|5] for a review of the top-down approach, [6H71I8] for the bottom-up approach and 
[9] for an hybrid approach). AdS/QCD is a powerful approach to the strong coupling regime 
of QCD because the string models depend on very few parameters (the string length, string 
coupling and number of colors and flavours of the dual theory). On top of that, general 
properties of hadron phenomenology (like vector meson dominance) seem to be universal 
in this approach in the sense that they do not depend on the particular AdS/QCD model. 
There is one string model, though, that has a field content very similar to large QCD in 
the regime of large distances and massless quarks. This is the Sakai-Sugimoto model [TOpT| 
that realizes confinement and chiral symmetry breaking. 

In this paper we investigate the production of baryon resonances with positive parity in 
proton electromagnetic scattering within the Sakai-Sugimoto model. First we define the cur- 
rent matrix element that describe the transition from a proton to a positive parity baryonic 
resonance. We derive the general expansion of the current matrix element in terms of scalars 
that we define as generalized Dirac and Pauli form factors which include as a particular case 
the well known elastic Dirac and Pauli form factors. Using holographic techniques we obtain 
a relation between the generalized Dirac and Pauli form factors and the couplings between 
baryons and vector mesons. This result can be interpreted as a realization of vector meson 
dominance in electromagnetic scattering of baryons. 

We also estimate in this paper the contribution of resonance production to the proton 
structure functions. The latter are Lorentz invariant scalars defined in Deep Inelastic Sca- 
tering (DIS) which is the inclusive scattering of a proton by a virtual photon (emitted by a 
lepton). The proton structure functions have a nice interpretation in terms of parton distri- 
bution functions that describe probability densities of finding partons (valence quarks, gluons 
and quark- ant iquark pairs) with a fraction of the longitudinal momentum of the proton. A 
calculation of these quantities at low momentum transfer proves difficult in non-perturbative 
QCD and usually relies heavily on some input from experiments or simulations. It should be 
noted that in the Sakai-Sugimoto model we can only make reliable predictions about scat- 
tering processes with low momentum transfers ( in this paper < 5{GeV)'^). This regime 
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is very far from the Bjorken limit of DIS {q^ — t- oo) which is well described by perturbative 
QCD. Moreover, it is beyond the scope of this paper to study the full inclusive DIS process 
with arbitrary final states. Therefore we limit ourselves to the contribution coming from 
single final state baryons with the same spin and isospin as the proton but different masses. 

Electromagnetic form factors have been obtained previously using bottom-up (phenomeno- 
logical) models and top-down string models. The meson form factors have been calculated 
in p!^15lll6j (bottom- up) and [T7|18fl9|l20] (top-down). Baryon form factors have been ob- 
tained in [2T|22] (bottom- up) and [23|24f25ll26] (top-down). Deep Inelastic Scattering in 
AdS/QCD was first investigated by Polchinski and Strassler in a bottom-up model for the 
case of scalar glueballs and baryon-like fermions [27] . Further development of DIS in bottom- 
up and top-down models include the large x regime [ 28|29|30ll3T] as well as the small x regime 
where Pomeron exchange dominates [32|f33|l34 .35 3 6,1^ . DIS structure functions have also 
been calculated for strongly coupled plasmas [38|39|40|l^ . 

The discussion in section |2] is independent of a specific model realization and contains 
novel results that apply very generally to non-elastic electromagnetic scattering for baryons, 
namely a detailed derivation of electromagnetic current matrix elements and structure func- 
tions. In section [3] we briefly review the Sakai-Sugimoto model and the description of holo- 
graphic baryons in this context. Section H] contains the derivation of the electromagnetic 
current and the generalized Dirac and Pauli form factors. Numerical results for the wave 
functions, masses, couplings and form factors are presented for the lowest excited states of 
spin 1/2 and positive parity. In section [5] we present our numerical estimate on the proton 
structure functions and the associated Callan-Gross relation. 

2. Baryon resonances in proton electromagnetic scattering 

Massless QCD with A'^^ flavours enjoys chiral symmetry {[/{Nf)^ x U{Nf)ji) at high ener- 
gies. At low energies chiral symmetry is spontaneously broken and the residual symmetry 
is vectorial corresponding to the group U{Nf)y. In the case of two flavours, the vectorial 
current can be written as Jy^ where a = (0, 1, 2, 3). The electromagnetic current JT"^ can be 
obtained as a combination of the isoscalar current Jv'" and the isovector current J,/'^ : 



In this paper we are interested in the production of baryon resonances with positive parity in 
electromagnetic scattering of a proton. Namely, we study the transition of a proton (denoted 
by B) with spin 1/2, isospin 1/2 and momentum p to a baryonic resonance (denoted by Bx) 
with the same spin and isospin as the proton but momentum px and different mass. This 
transition is characterized by evaluating the matrix elements of electromagnetic currents 
between the baryonic initial and flnal states. We denote the spin and isospin projections of 
the initial (final) baryon state as s {sx) and J3 (I^) respectively. 




a=0 



(2.1) 
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2.1. The generalized Dirac and Pauli form factors 



Evaluating the isoscalar and isovector current operators in the baryonic states we obtain 
current matrix elements that can be decomposed as 



2(27r) 



where 



{px - vY 



1 



5 



Kb 



;i,2,3) 



(2.3) 



and 0"" are the Pauli matrices. In fl2.2p we have used a generalization of the Gordon decom- 
position identity 



UBx{p'h^UB{p) = ubx{p') 



p'/^ + p^' iaf"" {pI - p^) 



rriB + mB^ 



rriB + ms. 



Ub{p) 



(2.4) 



The scalars Fg^^[q^) .F^^^^q^) in (12. 2p are the Dirac and Pauli form factors while F^^^{q'^ 
is required by current conservation 

= q,{Jrm - n{px, sx) [{px - p),l'Fj^^^{q') + iq' F'i'^^{q')\ u{p, s) 



u{px-,sx) i{mBx — mB)F^'^^ 

"^1 F^^^{q )\u{p,s) 



so that 



BE' 



[ruBx -mB)F^'^^{q 



(2.5) 



(2.6) 



The current matrix element now takes a transverse form 



{px, Bx, sx| ^(^'"(0)1^, B, s) 



2{2tv) 



+ KBa^xq^Fj^'Siq') u{p,s) 



u{px,sx) IvF^sAfl 



(2.7) 



We have chosen the signature rj^y = diag(— , +, +, +). The relativistic Dirac spinors that 
represent the initial and final states can be written as 



u{p, s) 



fXsip) 



, u{px,sx) 



X 



fx Xsx (P 



Xj 



Px 



fx 



-XsxiPx] 



(2. 



where / = y/E + itib and fx = \/ Ex + ttlbx- The two-component spinors Xs{p) and 
XsxiPx) are defined as eigenstates of the helicity operators 



p ■ a 



Xs{p) = sxs{p) 



Px ■ <y 

\px\ 
4 



XsxiPx) = Sx XsxiPx) 



(2.9) 



The convention for gamma matrices is 
1 



7 



7 



a' 
-a^ 



-1 

The baryon states are normahzed according to their charges as 



(2.10) 



(px, Bx, sx, If m\P: B, s, h) = {px, Bx, sx, I Jr(0) b, B, s, h) {27ryS'{q) 

= \{T'')ix,^5\q)5s,,5BB^F^S{^)- (2.11) 

In the case where the initial state is a proton and the final state has the same isospin polar- 
ization {I^ = = 1/2) the isoscalar and isovector charges are Nc/2 and 1/2 respectively. 
The elastic Dirac form factors are fixed at = as 



so we find from (12. lip 

{px,Bx,sx,l/2\p,B,s,l/2) = 5^{px-P)5sxs5bxB- 
The following spinor relations are very useful 



u{px,sx)l^u{p, s) 



xLipx) 



X 

^ ( f 



ffx + -rrP^ ' ^P' ^ 
J Jx 



Xs 



IX 



r2 , SXS\PX\\P\ 

Jx p 



xLipx)xs{p) 



(2.12) 



(2.13) 



(2.14) 



u{px, sx)l'u{p, s) 



2^EE^ 
^ ( f 



2VEE^ \fx 



j^s\p\ + sx\px\ 



xL{Px)cr'Xs{p) 



(2.15) 



u{px,Sx)(T°'qMp,s 



2^/EE^ 

i 



J Jx 



2^EE^ \fx 



^s\p\ - sx\p 



fj 
p 



x\ 



xlx{Px)(Ti(TjXs{p) 

qixlApx)o-\s{p) , 



(2.16) 



u{px, sx)(j'^qju{p, s) 



1 



2y^\fxJ 



fxfcTk 



{p + qTp^daakCTb Xsip) 



fxf 



'2 sxs\px\\p\ 

ix p 



xLiPx)(^kXsip) ■ 



(2.17) 
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The Breit frame. The Breit frame is characterized by the condition Ex = E which means 

that the photon has zero energy (go = 0). As shown in the appendix, choosing the photon 
in the z axis, the the photon and baryon momenta in the Breit frame take the form 



P3 



(0,0, 0,g) , 
, E 

2x 



(E,0,0,P3) , p^x=P^ + <i' 



(2.18) 



Using eq. fl2.7p and the spinor relations f l2.14p - fl2.17p we can work out the components of 
the current matrix elements. In the Breit frame the current matrix elements simplify to 



{px,Bx,sxJ^\fv"mP.B,s,h) 



P.a 



(2.19) 



(j>X:Bx,sxJt\JTmP:B,s,h) 



^(^si^")i^i-/'''(ljXlAPxHXs{p) 



X 



P,a I 2\ 



where 



a 

fi 
E 



^JE- 


f TUB 




- rriBx 


VE- 


f rriB 




-niBx 


g2 





[E + + {E- mB){l - 2x)] 

E + ITiBx 



2x 



E + ms 



+ 2x- 1 



To get fl2.19p and (12. 20 p . the following identity proved useful, 

sxlxi.Px)(y\s{.p) = -S'^xlx{px)xs{p) - i^^^xlx{.Px)(y^Xs{.p) , 

which is valid in the Breit frame only. 

Note that in the elastic case mBx = f^B, fx = f and kb = 1/(2^^5) so we obtain 

X\ jG,a, 



{px,Bx,sxjnJTW\P^B,s,h) 



(2.20) 



(2.21) 



(2.22) 



^ 2(^^"")^3-^3 (^) xUpx)Xs{plGr{<l') 

{px.Bx.Sxji\JTmP:B,sJ,) = (^-^y''%xUPxWxsmf^{q'^ 

(2.23) 



where 



D,a/ 2\ 



'—F'B\<f) 



Am 



G 



F^^\q') + F^^\q 



(2.24) 



B 



are the elastic electric and magnetic form factors also known as the Sachs form factors. 
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Figure 1: Exemplary diagram for a deep inelastic scatter- 
ing process. A lepton £ exchanges a virtual photon with a 
hadron of momentum p. 



2.2. Deep Inelastic Scattering and the proton structure functions 

Deep inelastic scattering (DIS) is a primary tool to investigate the internal structure of 
hadrons. DIS refers to the scattering process of a lepton on a hadron. The lepton interacts 
with a hadron of momentum pf^ via a virtual photon of momentum (cf. figure 1). The 
final hadronic state is denoted by X and momentum p'^. Such a process is commonly 

2 

parametrized by two dynamical variables, namely the Bjorken parameter x = and 
the photon virtuality (for a review of DIS, see e.g., [Hj). The standard limit in DIS 
corresponds to the Bjorken limit of large and fixed x. In this paper we are interested in 
the regime of small where non-perturbative contributions are relevant. 
The DIS differential cross section is determined by the hadronic tensor, 

= / d'xe'''-%p,s\ [j^{x),r{S))\ b,.), (2.25) 

where J^{x) is the electromagnetic current. Inserting the sum of the final states X we can 
write the hadronic tensor as 

^"' = ^y. y^C^^Y^^P + q - Px){p, s\J^mX){X\J''{0)\p, s) . (2.26) 

s X 

The structure functions Fi(x, g^) and F2{x,q'^) are Lorentz invariant scalars that appear in 
the decomposition of the hadronic tensor : 



W^-' = F,ix, q') (v"" - ^ ) + ^^2(x, q') + £ ) (^^ + £ ) • (2-27) 



2.2.1. The contribution from resonance production 

First of all we need to transform the spinors and baryon states of the previous section as 

uip.s) —^u{p,s) , \p,B,s) , ^ -|», (2.28) 

1/^' ' / y2E(27r)3/2 ' " ' 



in order to get the standard relativistic normalizations 



u{p, s)u{p, s) = 2mB , {px, Bx, Sx\p, B, s) = 2^ EEx{2Ti f5\px - M 



Using ([21]), (ETD and fl2:28|) we obtain for Jg 



{px,Bx,sx\J''mp,B,s) 



rx 



1/2, 



I 7]' 



u{px,sx) IvFbbM 



where 



(2.29) 



(2.30) 



are the Dirac and Pauh electromagnetic form factors. 

The baryonic tensor for a spin 1/2 baryon in the case where one particle is produced in the 
final state can be written as 



-yy 



s,sx ms. 



d^px_ 
(27r)3 



X {27r)^6\p + q-px){p,B, s\J^\0)\px, Bx^sx) {px, Bx, sxU^mp, B) 
= \j2J2^[(P + 9y + <x] (P' B, slJ'^iO) \px, Bx, sx) {px, Bx, sxU'iO) \p, 5, s) . 

s,sx ruB^ 

(2.32) 

Substituting fl230D into (12:321) we obtain 



Y^5[{p + qY + ml^] 



q^q" 



FEbx {<i')FEbx + F^Bx i<l')FEBx 



2n i?P 



+ F^Bx {<f)FEBx iQ")Cp. + Fl^Bx i<l')F^Bx il'^n 



pa 



(2.33) 



where 



A 



per 



-p^ip + qY tr(7^7p7^7<^) + ruBniBx tr(7p7<x) 

4| [mBrUBx + P-{P + q)] Vpo- - '^PpPa - PpQa - PaQp^ 



2 T 

mBiTiBx + P - TT 
2x 



Vpa - 2PpPa - PpQa - PaQp^ , (2.34) 



per 



-KBq q [mBniBx tr((TApfTA<,) - p^(p + qf tr(aAp7rf7Aa7r)] 
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-(qpPa + qaPp) + 
X 



rriBmB^ 



SI 

2x 



X 

^ -Yx 



g2 



(2.35) 



C 



pa 



and 



-zkbQ^ [rriBip + qY tr(crAp7r7a) + "i^^-p^ ii^axplalr)] 

4kb| [-mB{p + q)-q + niB^p ■ q] r]pa + ["^b(p + q)p - ^B^Pa] g^} 



4:KBq" 



Vp> 



qa ~1 

,^ + [niBqp + {rriB - mBx)Pp] ^|,(2.36) 



V 



pa 



iKBq^ [mB{p + qY tr{axalplT) + mB^P'' ^^{(^Xalrlp)] 

c 

^ap 1 



and we used the sum over spin formula 

u{p, s)u{p, s) = -il^Pp + niB 



(2.37) 



(2.38) 



and the gamma trace identities 



tr(7M7;^) 

^'^irip.lulpla) 

tr(7/.7!^7p7a7A7r) 

tr(a^^7p7a) 
iT{af,u(ypa) 

tr(cr^i.7pO-<^AO-r) 



+K'npvVpa - V^ipVva ~l~ ^pa^up) — ^Vpi'pa i 

'QpuVpaXr ^pp^vaXr ~\~ ^pa^vpXr ^pX^upar ~\~ ^jpr^vpaX 

- tr(crpi.7^7p) = M{-r]ppr]y„ + r]p^r]^p) , 
-4(-?7pp?7^a + VpaVi^p) : 

-4 - r]pp{r]^^riXr - rji^xrjar) + VpaVupXr - rjpxfjupaT 

-^r]pr{-r]uar]px + ?7i/A?7pa) 



(2.39) 



Note that the terms with qp or will vanish when contracting with the transverse tensors. 
Using flCTj) . (12:351) ■ (12361) . (fa?!) in (12331) we obtain the baryonic tensor 



pv 



q^q" 
q2 



Fi{q\x) 



2x 



(2.40) 



in terms of the structure functions 

Fi{q\x) = J] 5 [(p + g)' + m|^] 



-mBirriBx - "^b) + — 
2x 



+ 

+ 2 



q^ 



rriBiniBx + "2^) + — 1 

2x \x 



rriB + -^i^Bx - ^b) 



^lq\FEBx{q')f 
^Bq'FEBxiq')FEBxi^' 



(2.41) 
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and 



Interestingly, we can rewrite Filq"^, clS db binomial squared, i.e.. 



(2.42) 



(2.43) 



where 

C = q {niB + thb^ 



-1/2 



rriB + —{mBx - mB) 
2x 



1/2 



[{mBx - rriB) + q \ 



(2.44) 



Note that in the elastic case ( = q/V2 and kb = l/{2mB) so that the structure functions 
reduce to 



Fiiq^x) 



6{q' + 2p.q) [FEBiq') + F^Bm'- 



F2{q\x) = q^6{q^ + 2p-q) 



iFEsi^l')? + ^{FU<l')r 



Am 



B 



(2.45) 



2.2.2. The helicity amplitudes 



It is interesting to compare the result that we have obtained for the structure functions with 
the standard result in terms of the helicity amplitudes |15] 



Fl(g^x) = $^5[(p + g)^ + m|Jm|(G+^^(g^))^ 



X [{G^BBAq')r+2{G%BA<i')n 



4m^x^ 



(2.46) 



The helicity amplitudes G'^^^(g^) and G'^^^(g^) describe transitions when the initial state is 
a proton and the final state is a baryon of the same spin but different helicities. The amplitude 
G+ (Go) corresponds to a final baryon with helicity 1/2 (—1/2) and spin polarization —1/2 
(-1/2). 

Comparing (12.461) with our results for the structure functions we get the interesting relations 



c 



iG^BBA<i')? = \{fEbM) + fEbA^')? 



m 



{gIbM')? 



b 

q^\ C frriB^+mB D / 2 

-^BBx^^ 



2 J m% 



FEbA^")] (2-47) 



q^ 



niBx + rriB 
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3. Holographic baryons in the Sakai-Sugimoto model 

As pointed out in the introduction, the Sakai-Sugimoto model provides new insight into the 
problem of hadronic scattering in the non-perturbative regime. Moreover, baryons have been 
successfully incorporated into this model by several groups |^23f47|l2^ . This development 
was inspired by the Skyrme model [H] and Witten's original proposal of baryon vertices |49j . 
In this section we briefly review the Sakai-Sugimoto model and describe the construction of 
holographic baryons. 



3.1. Review of the model 

The Sakai-Sugimoto model is based on a configuration of D4 branes and A'^^ D8-D8 branes 
in the limit of large Nc with fixed Nf. This limit allows a supergravity description and can 
be interpreted as the quenching limit in QCD. The Sakai-Sugimoto model is the first string 
model that realizes confinement and chiral symmetry breaking. Below we describe this model 
in some detail. 

Consider a set of Nc coincident D4-branes with a compact spatial direction in type IIA 
supergravity [50]. The D4-branes generate a background with the following metric, dilaton 
and four- form: 

,,3/2 p3/2 J„,2 

f{u) = 1-- , e'^ = 9s^, , F4 = — ^64, (3.1) 

where is the tip of the cigar geometry generated by the D4 branes and R = {TrgsNcY^^y/a'. 
The r coordinate is compact and in order to avoid conical singularities the r period is fixed 
as 

47r i?3/2 

Sr^-^. (3.2, 

As a consequence, we get a 4-d effective mass scale 

M- — --^ (3 3) 

The r compactification is introduced as a mechanism of supersymmetry breaking and con- 
finement. Imposing anti-periodic conditions for the fermionic states we get at low energies 
a four- dimensional non-supersymmetric strongly coupled U{N^ theory at large with 't 
Hooft constant given by 

A = gl,,Nc = (27rM,) QsN, h . (3.4) 
It is convenient to define a pair of dimensionless coordinates y and z defined by the relations 



u = u. 



^{l+y^ + z') ^ u^Kli! , r = ^ arctan (^-^ . (3.5) 
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In terms of these coordinates the metric takes the form 



5/6 



(3.6) 



Now consider Nj coincident D8-D8 probe branes hving in the background generated by 
the Nc D4-branes. The probe approximation is guaranteed by the condition Nf <^ N^.. The 
Nf D8 branes introduce quark degrees of freedom as fundamental strings extending from the 
D4 branes to the D8 branes. The dynamics of the D8 and D8 branes is dictated by the DBI 
action. It turns out that the solution to the DBI equations smoothly merges the D8 and D8 
branes in the infrared region (small u). This is a geometrical realization of chiral symmetry 
breaking U{Nf) x U{Nf) — U{Nf). In the simplest case the solution is just y = (antipodal 
solution) and the induced D8-D8 metric takes the form, 

dsDs = ul/'R-'/'Kl/' r,,AxW + ^R"\l'^Kf'^dz^ + R:"\1'^ Kl'^d^l , (3.7) 

where = I + z"^ . Considering small gauge field fluctuations A^.Az depending only in 
and z directions the action of the D8-D8 branes reduces to a five- dimensional U{Nf) Yang 
Mills-Chern Simons (YM-CS) action in a curved background. The action of the model reads 



S — Sym + "S'cs 



•S'ym = —k / di^xdz tr 



where k = AiVc/(2167r'^). Here, fi,^ = 0,1,2,3 are four- dimensional Lorentz indices, and 
z corresponds to the fifth dimension. The quantity A = Aadx" = A^dx^ + Azdz {a = 
0, 1, 2, 3, z) is the five- dimensional U{Nf) gauge field and J-" = \j^apdx" /\dx^ = dA + iAAA 
is its field strength. 

3.2. Vector meson dominance 

The gauge field A^{x, z) can be expanded, in the Az = Q gauge, as 

oo 

Af,{x, Z) = %{X) + A^{x)ipo{z) + [v';^{x)ip2n-l{z) + a''^{x)ij2n{z)] , (3.9) 

n=l 

where 

V.(^) = lu~'[Ai + d,]U + \u[A^ + d,]U^\ 
M^) = lu~'[Ai + d,]U-^U[A^ + d,]U~\ 
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U{x) = , Al^''\x) = Alix) ± A^^{x) , 



(3.10) 



and the 'ipn{z) modes satisfy 

K j dzK~^'''^n{z)^m{z) = Km , 'K'^'d, \K,d,^^{z)\ = X^M^) ■ (S-H) 

Using the Kaluza-Klein expansion fl3.9p and integrating the z coordinate we get a four- 
dimensional effective lagrangian of mesons and external U{1) fields. The vector (axial vector) 
mesons are represented by the fields v^{x) {ci^{x)) and correspond to the modes ip2n-i{z) 
i4'2n{z)). The pion is represented by the field 7r(x) and corresponds to the mode ipoiz). In 
addition, we have external U{1) vector (axial) fields represented by A^ (A^). 

In order to have a diagonal kinetic term, the vector mesons are redefined as v"^ = + 
{gv"/M^n)V^ and the quadratic terms in the vector sector take the form [TT] : 

n 

where 

Ml„ = \2n^iMl , (7,n = K Ml„ j dz K:^l^i,2n-l{z) . 

The mixed term Qv^v^V^ represents the decay of the photon into vector mesons which is a 
holographic realization of vector meson dominance. 



Tr 



-V. 



3.3. Holographic baryons 

We restrict ourselves to the case Nf = 2. The U{2) gauge field A can be decomposed as 

A = A + A^ = A^^ + A^ = J2^''Y^ (3.12) 

a=0 

where r* (i = 1, 2, 3) are Pauli matrices and r° = I2 is a unit matrix of dimension 2. Thus, 
the equations of motion are given by 



647r2' 



-kK d F^" -\ ^^!f_f2A'2.../i5 (pa pa , p p ^=n 

-kK fV F'^T + g^A'2...M5 pa p =Q 

where V a = da + iA^ is the covariant derivative. 



13 



The baryon in this model corresponds to a sohton with a nontrivial instanton number. 
This is explained very well in the original paper by Sakai and Sugimoto [10], cf. also [16]. The 
reason for this is the following: The Sakai-Sugimoto model is a holographic model of large 
Nc QCD and naturally incorporates the Skyrme model. It has been known for quite some 
time [12] that the Skyrme model provides a qualitatively correct effective theory of baryons 
in large QCD. From symmetry considerations we can approximate the skyrmion solution 
by an instanton solution [13] , due to the fact that the symmetry group of the Skyrme model 
is SU{2) and a skyrmion corresponds to a topological soliton solution of the pion field with 
homotopy 713(6'^). Thus, the holographic dual of a skyrmion can be approximated by an 
instanton, i.e., a self-dual solution of the Yang-Mills equations in Euclidean space, which in 
our case is the four- dimensional space parameterized by (M = 1, 2, 3, z). 

The instanton number is interpreted as the baryon number Nb, where 

= j d^xdz eMiM2M3MAFlhM2^M3Mi ■ (3.13) 

The equations of motion are complicated nonlinear differential equations in a curved space- 
time, so it is difficult to find a general analytic solution corresponding to the baryons. 



3.3.1. Classical solution 



Since we are working in the large A regime, we can employ a 1/A expansion. It can be 
easily observed that Scs will be sub leading compared to S'ym, and therefore the leading 
contribution to the instanton mass comes from the YM action. As it turns out [lEj, it is 
possible to focus on a small region around the center of the instanton at 2; = (because 
the instanton size will scale as A~^/^), where the warp factor is approximately 1. The 
corresponding field equations will be solved by a BPST instanton with infinitesimal size 
p — )■ 0. As explained in ^6], including the contributions to the field equations from the 
CS term will induce a non- vanishing f/(l) electric field Aq and will stabilize the size of the 
instanton (determined by the minimum of the effective potential for p) at a finite value. The 
classical solution near 2; = corresponds to a static baryon configuration and is given by 



A 



M 



Nc 1 



(p2 + 



Aa = A 



M 



0. 



(3.14) 



with the definitions 



{z- Z)-i{x-X)-r 
I 



i 



zy 



x\ 



(3.15) 



where X^^ = {X^ , X'^ , X^ , Z) = {X,Z) determines the position in the spatial direction. 
The effective potential for p and Z can be calculated by taking into account the nontrivial 
^-dependence of the background (through K^) at order A~^ and reads 



p2 Ar2 1 
Mp,^)-Mo(1 + ^ + ^^ 



Z^ 



(3.16) 
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where Mq = Stt'^kM^ is the minimal (groundstate) mass of the baryons. The effective 
potential is minimized at 

3.3.2. Quantization 

The quantization of the solitons is facilitated by employing the moduli space approximation 
method to study a quantum mechanical problem on the instanton moduli space. For a more 
detailled discussion, the interested reader is referred to refs. Here we merely present 

the results for the wavefunctions and energies of the lowest excited baryon states in the 
slowly moving (pseudo-) moduli approximation. These (pseudo-) moduli are: 

X\t) , Z{t) , pit) , a\t) , (3.18) 

where X* and Z represent the center-of-mass position of the soliton, while p is the size of 
the instanton and [I = 1, 2, 3, 4) determines the orientation of the instanton in the SU (2) 
group space, with the condition {a^Y = 1. The SU{2) gauge field takes the form 

Am = VAt,V-' - iVOmV^' , (3.19) 

where is given by eq. f l3.14p and V satisfies the Gauss law constraint 

- tv-'v = -x^< + x^'fiOgYd'' > (3-20) 

with 

X"" = -i tr(r''a-^a) , a = 04 + iaaT" . (3.21) 

Inserting (13.191) into the effective action we get the Lagrangian of collective motion of the 
instanton 

, M, ^ 



° {X' + Z') + Mo J2 (pa/) (pa,) - Kff(p, Z) . (3.22) 
1=1 



Quantizing this system we find the baryon wavefunctions as eigenstates of the Hamiltonian. 
The relevant quantum numbers are B = {I, l3,np,nz) and its spin s. For example, baryon 
wavefunctions with Bn = (1, +1/2, 0, n) can be written as 

\Bnt)(xR{p)i/jBAZ){ai + ta2), (3.23) 

where 
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Note that eq. fl3.23p gives a representation of a baryonic state in terms of the collective 
coordinates fl3.18p of the moduli space of the instanton. The baryon masses can be easily 
gleaned from the relevant Hamiltonians and the resulting mass formula reads, 



M^M..,<i±i)!.A.,,H(l^^.M..?|. (3.5) 

One observation is in order: Since the photon couples to the baryons via vector mesons as 
a direct consequence of vector meson dominance and since the vector meson wavefunctions 
only depend on the coordinate z, it is clear that the intial and final state baryons must have 
the same p quantum number due to ort honor mality, while they may differ in the z quantum 
number, due to the additional contribution from the vector mesons to the relevant coupling 
constants etc. 



3. 3. 3. Extension of the soliton solution to large z 

The classical solution fl3.14p is valid only near z = Q. This solution can be extended to large 
z as long as we require p C, which is the condition of small size for the skyrmion. Under 
this condition the equations of motion linearize and the solutions can be found by defining 
Green's functions corresponding to the curved space generated by : 



G{x,z,X,Z) = nJ2M^)MZ)Yn{\x - X\) 



n=l 

oo 



H{x,z,X,Z) = nJ2^r.{z)MZ)Yn{\x-X\), 

n=0 



(3.26) 



where ipn{z) is the complete set of vector meson eigenfunctions, and (f)n{z) is another set 
defined by 



and Yn{r) is the Yukawa potential 



(f>n(,z) 



1 



-Aipn{z) (n=l,2, ...) 



Yjr) 



1 e-^"" 



An r 

The gauge field solutions found in [2l] for the case p <^ ^ can be written as 



A, 



~G{x,z,X,Z) 



2k 

A - ^ 

2n 



X 



2 

2 



■ d ■ d 
dXi dZ 



p d 



^f^ fx" d ^ p d 



2 V 2 dX'' pdZ 

2 2 a„-l 



A^ = 27rV{2iaa~^ + 27rVar''a 



pdX\ 
H{x,z,X,Z), 

X' ( e^^^ i -5^ 



G{x, z, X, Z) 



dXi 



dZ 



d 



dX" 



(3.27) 



(3.28) 



G{x,z,X,Z) 
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= -27i^p^Bir-Si-'^H{x,z,X,Z), (3.29) 
where 

A^ = AA^A-^ - iAd^A , A = ag-^V'K (3.30) 

4. Generalized baryon form factors in the Sakai-Sugimoto model 

Now we are going to extract the generahzed Dirac and Pauh form factors by comparing the 
matrix element of the vectorial current in (12. 7p with the one that can be obtained from the 
Sakai-Sugimoto model. The latter, denoted here by Jy^gs)^ '^^^ gleaned from holography 
as [21] 

Jviss) = -^{ }^ [^^^i] + [^^^i] } ' (4-1) 

where J-"^!, is the field strength associated with the classical field (I3.29p . When comparing the 
vectorial current of (12. 7p with the one in the Sakai-Sugimoto model we will use the following 
prescription 

r]^{px, Bx, sxl Jr(0)b, B, s) = 7]^{px, Bx, sx| J(^(55)(0)b, B, s) (4.2) 

where t]^ = {rjQ^ff) is the polarization of the photon and we choose to work with transverse 
photons satisfying the relation rj^q^ = in order to avoid the discussion of current anomalies. 

4.1. Electromagnetic currents in the Sakai-Sugimoto model 

Let us now represent the electromagnetic currents as operators in the moduli space of the 
collective coordinates of the instanton. Using (I3.29p and (14. ip one gets 



70,0 / N _ ^f-^ 

■fy[ss)i^) = 27r'K{pHT[T'^doi^T^S.-')]daHy + ^Gy 

- tT[r''aT^a-']X' [{dad, - 5,ad^)Hy + e'^'^djGy] } , 

fy\ss)i^) = -'27r\p^ti[T'aTaar^][id'd'' -S'''d^)Hy + e''''djGy] , (4.3) 



where 



Gy = -'^gyn^2n-l{Z)Y2n^l{\x - X\) , 

n 

Hy = -J2Y^dzi^2n^i{Z)Y2n^,{\x-X\) 



n 
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Mo ' 



i d _ P' 



and 



are the spin and isospin operators. Note that 



Z{d'Hv) - X'G 



V 



1 

Mo 



{d'Hv)Pz - GvP' 



Here we used the relation dzHy = —Gy- 
Defining the Fourier transform as 



'^v\ss)^^) ~ I ^^^^ ^ '^v\ss)i^) 



and using the identity 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



fe-*^F2n-l(|x-X|) 



-ik-X 



k"^ + X2n-1 



we find 



(4.8) 



(4.9) 



2 ' + \2n-1 



pi i 

Mo ^ Te^''"'^^'^" 



gv"dzip2n-iiZ) 



—d 

A2n-l(A:2 + A2„-l) [Mo ^ 



, (4.10) 



0,c 

y(55) 



gv"dz^p2n-l{Z) 



^ V + A2n-1 



A2n-i(A;2 + A 



_[! ^e^^'^PiA;,p2tr(r=araa-i) 

1 



i/tiP^ tr[r^9o(ar'a-^)] + — (P ■ A;^^ - PP^)p2 tr[r"ar'a-i] 

Mo 



(4.11) 



y(55) 



27r fee" 



-ik-X 



^Y^9v^^3l^-l{Z)_^ijaj^, 



P + A2n-1 



^ ^^na^^2n^i(Z) 

V A2n-l(^2 + A2„-l) 



2 jriaN 



p2tr(r^araa-i). (4.12) 



18 



Note that one term arising from Z cancels with another from X* and we have used the 
relation d\i\}n[Z^ ^ —XnipniZ). Now we calculate the expectation values of the Sakai- 
Sugimoto currents : 



{px, Bx, sx\Jv\ss)i^)\P^ ^' ^) 
We define the baryon states as 

\p, B, s, I3) = 
\px,Bx,sxJ^) = 



(2^ 



{px,Bx,sx\Ji}';;ss)0)\p,B,s). (4.13) 



1 



(27r)3/2 

1 

(27r)3/2 



e'P^-''\nB)K)\s,h) 



R ■ 



e^^^--'^M\n,)\sx,I^) 



R ■ 



(4.14) 



Here we make use of the results and definitions of a recent publication j^I], in which a 
relativistic generalization of baryon states and wavefunctions was discussed in detail. In 
particular, the spin and isospin part was defined as 



R 



1 f flsJs) 

f 



1 



3 \R 



X 



fx{sx,I^\ 



^x\Px[, jX 

7 \SX,h 

JX 



(4.15) 



where |s,/3) and {sx,Is\ are the non-relativistic initial and final states associated with the 
spin and isospin operators. Evaluating the currents in these states we get 



iJ 



0,0 

^(55") 



(0)) = j^,Y^sx,I^\s,h)nFhBAf)^ 



(4.16) 



i,0 



V{SS) 



m 



1 iVc 

(27r)3 2 



{sxJ: 



Mo Wn'^K " 



if) - Y^Fl^Af)W -<fnSa\\s,h)R, (4.17) 



(27r) 



{n,\{sx,inR{F'^BAf) 



^ + ^e->.g,p^tr(r^ar.a- 



FLAf) -w'tr[r^9o(ar^a-i)] 



R ■ 



(4.18) 



V{SS) 



m 



27r\ 



3, jria\ 



(27r)3 
X {np\p^\np){sxJ3 



me''\ + F^^Jq'){qY - fS 



R 



tr(r ar„a 



\S,I3) 



R ■ 



(4.19) 
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where 





n 




= E 

n 


Fls^if) 


= E 



{nBx\'4)2n~l{Z)\nB) 

+ A2n-1 
(?T-B^|<9z^2n-l(2')|nij) 



(4.20) 



A2„,-l(g^ + A2n-l) 

the momentum q is the photon momentum defined q = px — f> and we have used 

{px\e-"^-'^\p) = 5\k-p + px) . (4.21) 
The following relations are very useful and are presented here for completeness, 



{sxJ^\s,h)R 

(sx , I tr (r^ar„" ^ a) I s , J3 ) R 

{sxji\RSa \s,h)R 



1 



[ffx-'-^^)6,.,xUpx)Xs{p) 



ffx 

ssx\p\\px\ 



if fx - ^™^) r;.,3 xUpx) cr^Xsip) 



1 



if fx - '-^^) iT^)r-ixUpx) Xsip) , 



(//x-^^^)%,xL(pl)a,x 



Positive parity resonances in the Breit frame. In the Breit frame we get for positive 
parity resonances 

li^iifixlAp'x)Xsip)FlBxi^)^ 



(4'(5S)(0)) 

(4;55)(0)) 



2(27r)3 
2(27r)3Mo 



f3 

^ BBx 



1 

2^' 



- -,ae^^''q,aaF'BBM)\Xsi'P) 



(4(55)(0)) 



2(27r)3 



ir'')r-ixUPx)Xsi'P)FhBA^) 



Mo 



2(27r)= 



-Y ) (r^)/X/Klp>p)e^^%xL(Px)a„x.(p)i^iB,(g')(;4.22) 



where 



^ ' 2E 



JE + mB^ 



Mo 
3 



1+2W1+ 



\E + ms;, + (^ - mB)(2x - 1)] , 
_ 9i=\ 



AniB 



(4.23) 



and a is given in fl2.2ip . We have also used the relation 



(//. 



X 



ssx\p\\Px\ ' 
ffx ' 



/ 

— + - ssxiE - mB)\2x - 1|] 
Jx 



(4.24) 



and the identity fl2.22p which are valid only in the Breit frame. 
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4.2. The generalized Dirac and Pauli form factors 



Using the holographic prescription (I4.2p we can compare, in the case of positive parity 
baryons, the current matrix element of f l2.19p .( l2.20p with the Sakai-Sugimoto current matrix 
element in f l4.22p . As a consequence we get our main result: the Dirac and Pauli form factors 



r 1 



BB. 



C,a + Pa— 



2 

4Ayo 



a2 + p2q2 



(4.25) 



BBx 



1 

Kb 



Pi 



4Mo 



a2 + /32g2 



BB^, 



(4.26) 



BB^ 



ja + Paq' (f^) (p^) 

a2 + /32g2 



F'sbM') 



(4.27) 



1 

Kb 



/32g2 



(4.28) 



where a and /3 are given in fl2.2ip and ^ is given in fl4.23p . It is important to remark that 
the generalized Dirac and Pauli form factors depend on the collective coordinates of the 
instanton through eq. (14.201) and the expected value (p^). 
The electromagnetic Dirac and Pauli form factors read 



FEbM) 



BBx 



N, 



F'b'bM) + F^^Uq 



P,3 



(4.29) 



4-.2.1. The non-relativistic (large X) limit 
First of all note that 



TUBx = mB + -j^Bx^* , 



so that 



(4.30) 



X 



(4.31) 



Elastic case. In the elastic case we have ubx = 0, m^^ = ttib and x = 1 so in the large A 
limit we find 



KB V 4Mo 



2Mn 
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1 

Kb 



/3 



gi=i 



i + o 



(4.32) 



so that 



BB^ 



BBy 



l + O 



9i=o 



BBy 



l + O 



l + O 



(4.33) 
(4.34) 
(4.35) 



BBx 



9i=i 



l + O 



^BBx 



(4.36) 



Non-elastic case. In the non-elastic case, we have = 2, 4, 6, . . ., so that in the large A 
limit we have 



+ f3a 



rriB 
E 
_l_ 

2 



-1/2 



2n|^ 
3 g2 



KB V 4Mo 



Kb V 3 / 



2Afo ^ 

gi=Q TUB 



O 



XNr 



2 

TUB 



E V A 



1 + C 



AiV, 



Thus, in the non-elastic case, we obtain in the large A limit 



7^,0 (^2 
BBx 



XNr 



2 \ -1/2 



NcFhBAq') 



o{i) 



(4.37) 



(4.38) 



^P,0 (2 



BB. 



9i=o rnB 
2 E \XNr 



NcFhBxi<l') 



(4.39) 



7.1>,3 /„2 
BB. 



rriB 
E 



O 



FhBxi<l") 



(4.40) 



^^^,3 ,„2 



BBx 



91=1 



l + O 



XNr 



FhBxi^') 



(4.41) 
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4.3. Numerical results 



4-3.1. Baryon wavefunctions 

Here, we present the results for some low-lying baryon wavefunctions in fig. O With our 
framework, we restrict our presentation to the parity even baryon wavefunctions n = 2j, 
because, as we shall see below, the coupling constants under consideration involving an even 
baryon state (e.g. the proton) and an odd baryon state will yield zero. Furthermore, we 
are only considering baryon states with (np)initiai = ('^p)finai, which is zero in the case of the 
proton; all other possibilities will also produce vanishing results in the calculation of the form 
factors below. Tabled] summarizes the mass spectrum for the proton B = (1, +1/2, 0, 0) and 




Figure 2: (Normalized) wave functions "ifBiki^) ^^^^ even baryon states. 

its excitations Bx = (1, +1/2, 0, ra^). Note that we have chosen the proton mass = 
Mq = 940 MeV for obvious phenomenological reasons, although strictly speaking the masses 
are all proportional to Nc in the holographic limit Nc, A — ooll] 

^As noted in changing to a value of 500 MeV would result in a fairly realistic mass spectrum 
for the excited baryon states, cf. their table (5.35). However, changing M^, would alter the baryon wave 
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n 





1 


2 


3 


4 


5 


6 


7 


8 


mB„/GeV 


0.940 


1.715 


2.490 


3.265 


4.039 


4.814 


5.589 


6.655 


7.472 



Table 1: Some numerical values for the masses of the excited baryon states connected to the 
proton. 

4-3.2. Baryon form factors 

According to figure 1, we need to study the interaction of the baryons with (external) pho- 
tons. This process is represented by the electromagnetic form factors. The Dirac and Pauli 
form factors were discussed at length in section 2. Here we present our numerical results for 
the generalized baryon form factors. The infinite sums over vector meson states appearing 
in the mathematical description of form factors were approximated by including the first 
48 vector mesons states in the numerical computations. The wavefunctions ip2k-i{z) of the 
vector mesons were discussed at length in ref. [TU] . 

Some numerical results for these quantities are listed in table |2l It is interesting to compare 



k 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Mi 


0.6693 


2.874 


6.591 


11.80 


18.49 


26.67 


36.34 


47.49 


60.14 


y/^M2 


2.109 


9.108 


20.80 


37.15 


58.17 


83.83 


114.2 


149.1 


188.7 




5.767 


-2.610 


0.1902 


0.7664 


-0.5162 


-0.01955 


0.2118 


-0.08413 


-0.05348 


gv*'BoB2 


-0.9276 


2.4670 


-2.6239 


1.0560 


0.6404 


-0.9990 


0.2499 


0.3848 


-0.3049 


dv^'BoBi 


0.3655 


-1.2608 


2.0708 


-1.9409 


0.6966 


0.6716 


-0.9855 


0.2490 


0.4558 


Qv^'BoBq 


-0.1871 


0.7299 


-1.4596 


1.8556 


-1.3815 


0.1713 


0.8448 


-0.8371 


0.03738 


Qv^BoBfi 


0.1091 


-0.4595 


1.0352 


-1.5643 


1.5657 


-0.7918 


-0.3319 


0.9318 


-0.5634 



Table 2: Dimensionless squared masses and decay constants for vector mesons and coupling 
constants between vector mesons and baryons. 

our results with experimental data. For the vector mesons, the mass ratios for the first ex- 
cited states in our framework are m^ijmyx = 2.07 and m^s/rriyi = 3.14. These results can be 
compared with the experimental values [52] (m^2/m^i)exp = 1-88 and (7Ti„3/m„i)exp = 2.20. 
For higher excited resonances the discrepancy between our results and the experimental 
data increases. This should be expected since the spectrum of vector mesons follows a 
nearly quadratic Regge trajectory in this model [19]. The coupling constant describing 
the elastic interaction of a p meson with a proton can be compared with experimental 
data. In our framework we obtain Qv^BoBq = 5.77 that is consistent with experimental value 
Sv^'boBo ~ '^■^ ~ ^'^ '^^^ coupling constant that describes the decay of the resonance 

functions, vector meson decay constants, coupling constants and so forth. Therefore it is not permissible to 
merely adapt the mass spectrum by changing M*. 
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Pii(1710) into a p meson and a proton is gv^BoB2 = —0.93 in our framework. This value is 
compatible with the empirical value Is'^i^^^^l ~ obtained from the experimental decay 
width [54]. There is no available experimental data for higher order couplings. 

Using the results for masses and coupling constants we can now easily obtain the Dirac 
and Pauli form factors for the first few baryon states, as shown in fig. |3l Our results for the 
elastic form factors of the proton are compatible with the results obtained in chiral soliton 
models [55]. Again, this should be expected from the fact that baryons in the Sakai-Sugimoto 
model are holographic representations of four dimensional skyrmions. 

Observe that the form factors -^^0^4 negative definite. This is not problematic since 
they will only appear in bilinear combinations in the derivation of the structure functions 
below, so that the sign will not matter. 



4-3.3. Helicity amplitudes 



In figure H] we present our numerical results for the helicity amplitudes 

that have been discussed in section 12.2.21 As before, we study the h 
amplitudes in the large A limit, where the expressions simplify to 



\/2mB 



1/2 



(4.42) 



The same limitations as for the generalized form factors apply here as well, i.e., the numerical 
errors become significant for j > 4 and > 5(GeV)^. Therefore the increase observed in 
for j = 3,4 and greater than 4(GeV)^ may be an artefact of the numerics. 
Knowledge of the helicity amplitudes \G~^^Q^.{q'^)\ is sufficient to determine the helicity am- 
plitudes cf. e.g. [15]. The results are visualized in figure [51 The model at hand 

predicts A^^^^^,{q'^ — 0) = 0(GeV)~^/^, which is consistent with experimental data, cf. e.g., 
[521I3B] . The experimental result quoted in the PDG summary tables 



IS 



Ai/2(g2 = 0) = 0.009 ± 0.022(GeV) 



-1/2 



(4.43) 



for the transition between Pii(1710) and the ground state. Other theoretical, i.e., quark 
model predictions are available in the literature, e.g., by Capstick and Arndt et. al. [57|l58] . 



Ai/2{q^ = 0) = 0.013(GeV) 



-1/2 



(4.44) 



for the above mentioned transition. 
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Figure 3: Dirac and Pauli form factors FQ^_^Q^.{q^) for the first four baryon states. The 
elastic case corresponds to j = 0, while j > yields the transitional form factors. For 
j > 4, the numerical errors become relatively large; we can only trust our results for approx. 
<5(GeV)2. 
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Figure 4: Helicity amplitudes \G~^'^^^.(q'^)\ plotted versus in (GeV)^. The transitions from 
protons to the baryonic final states are labelled by j, where j = 1 (blue, solid), j = 2 (red, 
dashed), j = 3 (orange, dotted) and j = 4 (green, dotdashed). 
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Figure 5: Helicity amplitudes \A^^^Q^_,{q'^)\ plotted versus in (GeV)^. As before, the 
transitions from protons to the baryonic final states are labelled by j, where j = 1 (blue, 
solid), j = 2 (red, dashed), j = 3 (orange, dotted) and j = 4 (green, dotdashed). 



5. The proton structure functions 

5.1. Baryon Regge trajectories in the Sakai-Sugimoto model 

Assuming approximate continuity of the mass distribution, we can now approximate the 
delta distributions in the following way: 



rris^ - s\ 



dn 



dn 



-\ -1 



6{n — n) 
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dn 



-\ -1 



fin), 



(5.1) 



with the definition 



s:= -ip + qf = ml +q'^{--l). 



X 



(5.2) 



Therefore we have to evaluate the Regge trajectory of the baryon spectrum in order to 



calculate We find (cf. eq. (^M) 



dn 



(5.3) 



where Mq can be chosen to match, e.g. the proton mass m^o and n := nz- 



We obtain an asymptotically quadratic Regge behavior for the baryon resonances in ac- 
cordance with expectations from the quantized Skyrme model, which is a feature inherited 
by the Sakai-Sugimoto model, cf. eq. (3.5). However, this is different from the experimen- 
tally observed linear behavior of the Regge trajectory. This is a well-known shortcoming of 
many holographic models of vector mesons and baryons (see for instance [3)16] ). There exist 
holographic models that exhibit a more realistic linear Regge trajectory by considering a 
dilaton field in the bulk, e.g. the soft-wall models |^2E3] . 



5.2. Numerical results 

5.2.1. Dependence o/-Fi,2 on and x 

It is now possible to extract information about the structure functions Fi^2iq'^iX) employing 
the following strategy: From the discrete set of baryon mass states and the relation (15. 2p . 
we find 

^m%, ■■= - = 9^ " l) • (5-4) 

It is possible for each j > to extract a discrete set of values for for a given set of fixed 
Bj0rken parameters, e.g., x = 0.001, 0.01, 0.05, 0.1, 0.3. Alternatively, we may calculate a 
set of values for x for a given set of, e.g., = 0.1, 0.5, 1, 2, 3 (GeV)^. The corresponding 
values are summarized in table El Lastly, we need to collect some results about the calculation 
of the isoscalar and isovector magnetic moments for the states under consideration. For the 
proton and its excited states with rip = and spin up, we find (cf. eqs. (3.16) and (3.32) of 

/^ko = - 0.842 6'^ f^s, /iki = ^ "^^^lH^ Pel - 3.52 6^^ (5.5) 

measured in units of the Bohr magneton /i^ = l/{2mBa)- Here, the mass and size of the 
classical instanton are Mq = Stt'^kM^: = 558 MeV and p^j = g-s^^-^, respectively, and we 
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j 


f{n = 2j)/M2 


Am|^y(GeV)2 


gV(GeV)2 


X 


1 


0.233 


5.317 


0.00532 for x = 0.001 
0.0537 for x = 0.01 
0.280 for X = 0.05 
0.591 for X = 0.1 
2.279 for x = 0.3 


0.0185 for g2 = 0.1(GeV)2 
0.0860 for g2 = 0.5(GeV)2 
0.158for g2 = i(GeV)2 
0.273 for g2 = 2(GeV)2 
0.361 for g2 = 3(GeV)2 


2 


0.144 


15.430 


0.0154 for a; = 0.001 
0.156 for X = 0.01 
0.812 for X = 0.05 
1.714 for X = 0.1 
6.613 for X = 0.3 


0.00644 for = 0.1 (GeV)^ 
0.0314 for g2 = 0.5(GeV)2 
0.0609for g2 = i(GeVY 
0.115 for g2 = 2(GeV)2 
0.163 for g2 = 3(GeV)2 


3 


0.0814 


30.353 


0.0304 for X = 0.001 
0.307 for X = 0.01 
1.598 for X = 0.05 
3.373 forx = 0.1 
13.008 for X = 0.3 


0.00328 for q^ = 0.1(GeV)2 
0.0162 for g2 = 0.5(GeV)2 
0.0319for g2 = i(GeV)2 
0.0618 for g2 = 2(GeV)2 
0.0899 for q^ = 3(GeV)2 


4 


0.0436 


54.947 


0.0550 for X = 0.001 
0.555 for X = 0.01 
2.892 for x = 0.05 
6.105 for X = 0.1 
23.549 for x = 0.3 


0.00182 for g2 = 0.1(GeV)2 
0.00902 for q^ = 0.5(GeV)2 
0.0179for g2 = i(GeV)2 
0.0351 for g2 = 2(GeV)2 
0.0518 for g2 = 3(GeV)2 



Table 3: Some values for and x according to eq. (15. 4p . 



have set = 3 for obvious phenomenological reasons. Sometimes it will be more convenient 
to utilize magnetic gj factors, which can be defined as follows, 

= j^<^\ (5.6) 

where ex* are Pauli matrices. The numerical values in the Sakai-Sugimoto model turn out to 
be 

(77=0^ 1.684, (77=1^7.031. (5.7) 

It is now a fairly straightforward exercise to evaluate the structure functions -Fi^2(9^, a^) for 
the discrete set of values given in table [31 Again, it should be stressed that we work in the 
large A limit, i.e., we only keep the leading terms in the large A, large Nc expansion before 
setting the baryon masses (which are of order 0{\N^) to there phenomenological values. 
The results are presented in figures E] and [3 

5.2.2. Comparing our results with experimental data 

The results shown in Figures [6] and [7] indicate that our structure functions Fi 2(g^,x) seem 
to be at least two orders of magnitude smaller than what is expected from DIS experimental 
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Figure 6: Structure functions Fi^2{,'f') for x = 0.3 (orange, solid), x = 0.1 (red, dashed) and 
X = 0.05 (green, dotted). 



data (cf., e.g., [52], chapter 16). A possible factor that may contribute to the discrepancy 
between our results and experimental data is the naive approximation of the Delta distribu- 
tion that we used in eq. (15.11) . Alternatively, we can approximate the Delta distribution by 
a Lorentzian function [45] : 

-1 

(5.8) 

where F^^ is the decay width of the resonance i?x- Identifying the first positive parity 
resonance B2 with the experimentally observed Pii(1710) and using the empirical value 
= 100 Mev we obtain the results for the structure functions shown in Figure [HI Note 
that the structure functions have improved by an order of magnitude. Unfortunately, we 
can not follow this procedure for the higher resonances because there are no experimental 
results available for the decay widths. 

The results for the proton structure functions obtained herein are understood to be non- 
inclusive and only represent a small fraction of possible final states, namely single final state 
baryons (the excited states of the proton) with spin 1/2 and positive parity. If we include 
the contribution from final states with negative parity [59] as well as final states with higher 
spiiJl and pion production relevant in this kinematical regime, we should get a better picture 
of the proton structure functions. 



F^ 



5.2.3. Callan- Gross relation 

The Callan-Gross relation can be easily studied numerically in this framework in order to 
check its validity for spin 1/2 particles. We define the ratio 

Rcoiq^x) := /!,^fr\ > (5-9) 
2xFi(g^, X) 



^See [22] for transition to A resonances in the Sakai-Suginioto model 
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Figure 7: Structure functions -^1,2(2^) for = 3(GeV)^ (purple, solid), = 2(GeV)^ (blue, 
dotdashed), = l(GeV)2 (green, dotted), = 0.5(GeV)2 (red, dashed) and g^ = 0.1(GeV)2 
(orange, solid, barely visible) . 

and use the numerical results for the structure functions (Figures [6] and [7j). The result is 
shown in Figure [9] and indicates an approximate Callan-Gross relation, i.e. Rqg ~ 1 for 
intermediate values of the Bjorken variable x. 

6. Conclusions and Outlook 

In the present paper we have derived a generalization of the notions of baryon electromag- 
netic form factors to non-elastic scattering when baryonic resonances with positive parity 
are produced. Subsequently, we have used the Sakai-Sugimoto model to estimate these 
form factors in the non-perturbative regime of large Nc QCD. We have also estimated the 
contribution of these form factors to the proton structure functions. 

It is important to remark that, in this paper, we have considered the large A limit in the 
Sakai-Sugimoto model where the description of baryons as solitons is derived from classical 
small instanton solutions, whose size is of order ~ A~^/^. There may exist 1/A corrections 
that can be dominant at large distances as suggested in a recent analysis [60161] . One should 
also bear in mind the limitations of describing baryons in the large Nc limit □. Interestingly, 
in [63] holographic baryons were constructed for the (bottom-up) hard wall model as solitons 
in an AdS^ spacetime with cut-off. These baryons have finite size and have the expected long- 
distance properties [60]. It would be interesting to investigate baryon resonance production 
in this model and recent holographic models such as [M1I65] or the recent string models 
based on the singular [66)120] and deformed [671168] conifold backgrounds. It would be also 
interesting to estimate the contribution of resonance production in other scattering processes 
like dilepton production in proton-proton scattering (see [69]). 

■^There occurs a phase transition at iVc w 8 that separates the small Nc (= 3) regime where nuclear 
matter behaves like a quantum liquid from the large Nc (holographic) regime where nuclear matter behaves 
like a crystalline solid [SI] 
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Figure 8: Structure functions -^1^2(3^) for = 3(GeV)^ (purple, solid), = 2(GeV)^ (blue, 
dotdashed), = l(GeV)^ (green, dotted) near the first positive parity resonance B2 
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A. Breit frame for inelastic scattering 

Consider the scattering between a virtual photon and a hadron in the hadron rest frame. 
After two rotations we can set the spatial momentum of the photon to the direction, so 
that 

if = (m^j, 0,0,0) 

g'^ = (go, 0,0,^3), (A.l) 
and we choose > 0. The virtuality and Bjorken variable in this frame are given by 

2777,5 go 

Now we perform a boost in the direction so that 

p'^' = (77775, 0, 0, -/37mB) 

q''' = (790 - /St^s, 0, 0, -/37go + 7^3) • (A.3) 
The Breit frame is defined by the condition gQ = so that 
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Figure 9: Callan-Gross ratio Rcg{x) for = 3(GeV)^ (purple, solid), = 2(GeV)^ (blue, 
dotdashed) and = l(GeV)^ (green, dotted). 



and we arrive at 



with 



q'^ = (0,0,0,Q), (A.5) 



P3 = -f. (A.6) 
2x 
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